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Abstract 

We show that any non-degenerate vector field u in L°°(Q, R ), where fl is a bounded domain in WL N , 
can be written as 

u(x) = ViH(S(x),x) for a.e. (1) 

where S is a measure preserving point transformation on Q such that S 2 = I a.e (an involution), and 
H : M. N x — > R is a globally Lipschitz anti-symmetric convex-concave Hamiltonian. Moreover, it is 
a monotone map if and only if S can be taken to be the identity, which suggests that our result is a 
self-dual version of Brenier's polar decomposition for the vector field u as u{x) — X7(p(S(x)), where tp is 
convex and S is a measure preserving transformation. We also describe how our polar decomposition 
can be reformulated as a self-dual mass transport problem. 



1 Introduction 

Given a bounded domain SI in R N , a classical theorem of Rockafellar [9i yields that a single- valued map u 
from SI to R™ is cyclically monotone, i.e., for any finite number of points (xi)™ =0 in SI with x$ = x n , we have 

n 

J2(u{x k ),x k -x k -i) > 0, (2) 



if and only 

u{x) = Vip(x) on 17, (3) 

where <p : M. n -> R is a convex function. 

On the other hand, a result of E. Krauss [7] yields that u is a monotone map, that is for all (x, y) in Q, 

(x-y,u(x)-u(y)}>0, (4) 

if and only if 

u(x) = ViH(x, x) for all x € Ct, (5) 
where H is a convex- concave anti-symmetric Hamiltonian on R N x R N . 

More remarkable is the polar decomposition that Y. Brenier [3] establishes for a general non-degenerate vector 
field, and which follows from his celebrated mass transport theorem. Recall that a mapping u : fl — ¥ WL N is 
said to be non-degenerate if the inverse image 7i -1 (7V) of every zero-measure N C R N has also zero measure. 
Brenier then proved stating that any non-degenerate vector field u S L°°(f2,R Ar ) can be decomposed as 

u(x) = Vip ° S(x) a.e. in f2, (6) 
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with ip : M. — > R being a convex function and S : f2 — > fl is a measure preserving transformation. 

In this paper, we shall prove another decomposition for non-degenerate vector fields, in the same spirit of 

Brenier's, but which can be seen as the general version of Krause's. Indeed, here is the main result of this 

paper. 

Theorem 1.1 Let ft be an open bounded set in R N such that meas(dfl) = 0. 

1. If u G L°° (Q,M. N ) is a non-degenerate vector field, then there exists a measure preserving transfor- 
mation S : 0, — >• f2 such that S 2 = I (i.e., an involution), and a globally Lipschitz anti- symmetric 
convex-concave Hamiltonian H : R N x R N ->• R such that 

u(x) = ViH(Sx, x) a.e. x <E Q. (7) 

2. If u is differ entiable a.e., and the map 

x ->■ (u(x), y x - y 2 ) + (u(yi) - u(y 2 ), x) (8) 

has no critical points in f2 unless y\ — y 2 , then there exists a unique measure preserving involution S 
such that satisfies <\21\ for some convex- concave anti- symmetric Hamiltonian H . 

3. In particular, if u is a strictly monotone map, then S is necessarily equal to the idendity. 
Since S is an involution and H is anti-symmetric, one can deduce from (ll.l[) that 

u(Sx) = -V 2 H(Sx,x) a.e. x eQ,. (9) 

The fact that S is a measure preserving involution provides an improvement on the first factor in Brenier's 
decomposition. On the other hand, the second factor i.e., the potential Vip, is obviously better than the 
partial gradient of a convex-concave anti-symmetric Hamiltonian on R N x R N , which is only a monotone 
map. The connection to self-duality will be described later in this introduction. 
We now give a few examples of how this decomposition appears in concrete situations. 

1. Basic monotone operators not derived from a potential: If u = Vip + A, where ip is convex and 
A is a skew-adjoint matrix, then clearly 

H(x, y) = <p(x) - <p(y) - (At, y), (10) 

and S is the identity. This is of course a very important case of maximal monotone operators that we shall 
discuss later. 

2. Helmholtz Decomposition of vector fields: Let f2 be a smooth bounded connected open set in R N . 
Any smooth vector field u on D, can then be written in a unique way as u(x) — Vp(x) + v(x), where p is a 
smooth real function on f2, and v is a smooth divergence free vector field parallel to the boundary of f2. By 
considering u e as a smooth perturbation of the identity map: 

u e (x) = x + eu{x), x G fl. 

for e small enough, we can write 

u t (x) — ViH(x,x), for all x G Cl, 

where 

H(x,y) = -\x\ 2 + ep{x) + e(x - y, - ) - -\y\ 2 - ep(y), 

and again S is the identity operator. Note that for e small enough, H is convex in the first variable and 
concave in the second one. 
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Note that both examples above fit in the framework of the result of E. Krauss [7J, who -as mentioned above- 
has shown that if u is a single-valued maximal monotone map on M. N , then u(x) — ViH(x,x), that is S 
is the identity map whenever u is a monotone map. We shall come back to this theme when we discuss 
self-duality below. Now, we give examples of non-monotone operators. 

3. Decomposition of real matrices: Any N x N matrix A can be decomposed as A s + A a where A s is 
the symmetric part and A a is the anti-symmetric part. The symmetric part A s can then be written as the 
product RS of a symmetric non-negative matrix R and a real unitary matrix S. It is not difficult to check 
that since A s is symmetric so is S and therefore S 2 = I. 

It follows that we can write the following decomposition for the matrix A : 

Ax = A a x + RSx = ViH(S(x),x), for all x E R N , 

where 

1 1 

H{x,y) = -(Rx,x) - ~{Ry,y) - (A a x,y) 

is clearly a skew-adjoint Hamiltonian and S is a symmetric involution matrix. 

Note that the condition (J5|) is, in this case, equivalent to saying that the symmetric part A s of A is non- 
singular. Indeed, 

(u(x),yi - m) + (u(yi) - u(y 2 ), x) = 2{A s {y 1 - y 2 ),x) 

has a critical point if only if A s y\ — A s y 2 , which means that y\ = y 2 whenever A s is assumed to be non- 
singular. This is compatible with the classical fact that R and S in the above decomposition of A s are unique. 

4. Examples of representations on non-monotone maps on the line: 

(i) A simple non-monotone example is the function 

u(x) = sinx + xcosx. (11) 
It can be written as u{x) = ViH(S(x),x) on [0,7r], where 

H(x,y) = xsiny — yshix and S(x) — tt — x. (12) 

(ii) More generally, a large class of examples of convex-concave anti-symmetric Hamiltonians is given by 

H(x, y) = f(x)g(y) - f(y)g(x) + h(x - y), 

where h is odd and with suitable conditions that render H convex in x. For a G [0, 1], the function 

, , _ f f'(a-x)g(x) -g'(a-x)f(x) +h'(a-2x) if < x < a . . 

K > 1 f'(x)g(x)-g'{x)f(x) + h'(0) if a < x < 1. 1 1 

can then be written as u{x) — Vi-ff (S(x), x), where S(x) = a — x on [0, a] and S(x) — x on (a, 1). 

(iii) A more interesting example is the map 

, \ f 3 - 2x if i < x < 1 ,,,, 
U ^ = { 2x if 0ix< |. (U) 

One can easily verify that u(x) — ViH(Si(x),x), where Si(x) = 1 — x and H is given by the following 
formula 

!-2xy + 2x-y-\ if < x < \ and ± <y <l 

x — y if0<x,y<i or \<x,y<l (15) 

2xy - 2y + x + \ if | < x < 1 and < y < \. 

Also note that u(x) — \7iH(S 2 (x),x), where S 2 (x) = x + | on [0, i) and S 2 (x) = x — \ on (|, 1], which 

means that the involution S appearing in the decomposition is not necessarily unique. 

Actually, one has non-uniqueness whenever there exists two subsets A, B of positive measure such that 

u(x).y + u(y).x = f(x) + g(y) when (x, y) E A x B, 
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for some functions / and g. This is what happens in the previous example with A — [0, 1/2] and B = [1/2, 1]. 
It follows that x — > u{x).{y2 — yi) + {u(y2) — u(yi)).x doesn't depend on x. If u is diffcrcntiable then 
it' (a) + u'(y) = for (x, y) e Ax B. 

5. Why can the decomposition be considered "self-dual"?: Let X be a reflexive Banach space, and 
recall from [5] the notion of a vector field dL that is derived from a convex lower semi-continuous Lagrangian 
on phase space L : X x X* — > K U {+00} in the following way: for each x 6 X, the -possibly empty- subset 
dL(x) of X* is defined as 

BL{x) := {p 6 X*; (p, x) 6 0L(x,p)}. (16) 

Here dL is the subdifferential of the convex function L on X x X*, which should not be confused with 
91/. Of particular interest are those vector fields derived from self-dual Lagrangians, i.e., those convex lower 
semi-continuous Lagrangians L on X x X* that satisfy the following duality property: 

L*(p,x) = L(x,p) for all (x,p) €X X Jf*, (17) 

where here L* is the Legendre transform in both variables, i.e., 

L*{p,x) = sup{(y,p> + (x,q) - L{y,q) : (y,q) £ X x X*}. 

Such Lagrangians satisfy the following basic property: 

L(x,p) - {x,p} > for every (x,p) £ X x X* . (18) 

Moreover, 

L(x,p) — (x,p) — if and only if (p,x) £ dL(x,p), (19) 
which means that the associated vector field at x 6 X is simply 

dL(x) := {p e X*-L{x,p) - {x,p) = 0}. (20) 

These so-called selfdual vector fields are natural but far reaching extensions of sub-differentials of convex lower 
semi-continuous functions. Indeed, the most basic selfdual Lagrangians are of the form L(x,p) = tp(x)-\-(p*(p) 
where ip is such a function in X, and tp* is its Legendre conjugate on X*, in which case 8L(x) — dtp(x). 
More interesting examples of self-dual Lagrangians are of the form L(x,p) = ip(x) + Lp*{—Tx + p) where tp 
is a convex and lower semi-continuous function on X, and L : X — > X* is a skew adjoint operator. The 
corresponding selfdual vector field is then dL(x) = Tx + dip(x). 

Actually, it turned out that any maximal monotone operator A (a notion studied for example in J4 ) is a 
self-dual vector field and vice- versa [6]. That is, there exists a selfdual Lagrangian L such that A = dL. This 
fact was proved and reproved by several authors. See for example, R.S. Burachik and B. F. Svaiter [5], B. 
F. Svaiter 11 , and Baushke and Wang [?]). 

This result means that self-dual Lagrangians can be seen as the potentials of maximal monotone operators, 
in the same way as the Dirichlet integral is the potential of the Laplacian operator (and more generally as 
any convex lower semi-continuous energy is a potential for its own subdifferential). Check [6] to see how 
this characterization leads to variational formulations and resolutions of most equations involving monotone 
operators. 

Consider now the Hamiltonian Hl on X* x X* corresponding to L, that is the Legendre transform of L in 
the first variable, 

HL(p,q) = sup{(x,p) - L(x,q);x £ X}. 

It is convex-concave and satisfies Hi J (q,p) < —H^p.q). In most concrete examples, it is actually anti- 
symmetric. If now A is a maximal monotone operator, then A" 1 is also maximal monotone and therefore 
can be written as A^ 1 = dL, where L is a selfdual Lagrangian on X* x X that can be constructed in the 
following way: First, let 

N(p,x) = sup{(p,y) + (q,x - y); (y, q) £ Graph(A)} (21) 
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in such a way that 

N*(x,p) > N(p,x) > (x,p) for every (x,p) € X x X*. (22) 
Then consider the following Lagrangian on X* x X, 

L(x,p) := inf |^JV(px,a;i) + ^N*(x 2 ,p 2 ) + j^\\xi - x 2 \\ 2 + ^\\pi -P2W 2 ; (x,p) = ^(xi,pi) + ^(x 2 ,p 2 ) 

It was shown in [B] that L is a self-dual Lagrangian on X* x X. One can also show that the corresponding 
Hamiltonian Hl on X x X is anti-symmetric, and that 

p e Ax if and only if (p, — p) £ 8Hl(x, x). (23) 

Moreover, Ax = V \H{x, x) if ^4 is single-valued maximal monotone operator, which is Krause's result 
mentioned above. 

Compared to Brenier's, our decomposition now looks like we have replaced the potential of a convex function 
in Brenier's theorem by a more general maximal monotone operator A (or a self-dual Lagrangian L), while 
we have strengthened S to make it a measure preserving involution. 

6. Connection to Monge's mass transport: We shall see in the next section that the transformation S 
appearing in the decomposition (JJJ) of u maximises the quantity 

(u(x), S{x)) dx, 

n 

among all measure preserving involutions on S7. Equivalently, it does minimize 

\u(x) — S(x)\ 2 dx, 



Q 

which is the distance of u to the set of measure preserving involutions on SI. The latter minimization can 
now be related to an optimal transport problem with a quadratic cost, between the measure /x on SI x u(tt) 
obtained by pushing Lebesgue measure on fi by the map x —> (x,u(x)), and the measure v on u(fi) x fi 
obtained from p by the transposition map (x,y) — > (y, x). Indeed, any map T pushing p onto v can be 
parametrized by an application S : fi —> fi via the formula: 

T : (x,u(x)) — > (u(Sx), Sx), 

and the transport cost is then equal to 

- / [\u(Sx) - x\ 2 + \u(x) - S(x)\ 2 ]dx, 

which, in the case where S is a measure preserving involution, coincides with \u(x) — S(x)\ 2 dx. 
Note now that if T is an optimal transport mapping fi onto v, then the map (y, x) — > (x, y) — > T(x, y) would 
be an optimal transport mapping v onto p. It will then follow that if the optimal transport T from fi onto 
v was unique, then the S corresponding to T would necessarily be an involution. Now in terms of Brenier's 
theorem, the uniqueness would necessarily lead to a convex function L on x R N such that T = VL and 
where L is a selfdual Lagrangian, i.e., L(x,p) = L*(p,x). The anti-symmetric Hamiltonian H would simply 
be the Legendre transform of L with respect to the first variable. 

Unfortunately, the measures p and v on the product space are too degenerate to fall under the framework 
where we have uniqueness in Brenier's theorem. Hence the need to establish the result directly and without 
resorting to Mass transport. 

On the other hand, if one starts with a measure fi on the product space SI x SI that is absolutely continuous 
with respect to Lebesgue measure, then one can apply Brenier's theorem to find an optimal transport map 
VL that pushes fi onto its transpose p via the involution R{x,y) = (y,x), and where L is a convex function 
Consider now the convex function Mix, y) — L*(y, x) — L* o R(x, y) where L* is the Legendre 
transform of L with respect to both variable, and note that VM = R o VL* o R also maps p into p. By the 
uniqueness of the optimal map we deduce that VL = VM = Ro VL* o R, which means that L is a self-dual 
Lagrangian (up to a constant). 

We are extremely grateful to Bernard Maurey for very insightful discussions, and to Ivar Ekeland and Philip 
Loewen for their helpful input. 
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2 A variational formulation of the problem 



The proof of the standard polar factorization by Brenier was based on tools from the Monge-Kantorovich 
theory. Later W. Gangbo [5 provided a more direct proof by formulating the Brenier decomposition as 
the Euler-Lagrange equation corresponding to a suitable variational problem. Our approach is in line with 
Gangbo's method and involves various new results about skew-symmetric functions, which may have their 
own interest. 

To formulate our variational problem, we start by considering the set % of all continuous anti-symmetric 
functions on f2, that is 

H = {H e C(n x n);H(x,y) = -H(y,x)}. 

For each H £ H, define L H : fi x R N -> R by 

L H (x,p) = sup{(t/,p) - H(y,x)}, (24) 
yen 

where (., .) is the standard inner product in R N . Set C = {L H ;H £ %}. 
Here is our main theorem. 

Theorem 2.1 Let ft be an open bounded set in W N such that mea(dfl) — 0, and let u £ L°°(f2,K ) be a 
non- degenerate vector field. Consider the following two variational problems: 

P 00 =mt{j n L H (x,u(x))dx;H £ U) (25) 

and 

Doc = sup { J n (u(x), S(x)) dx; S : ti ^ is a measure preserving involution} . (26) 
Then the following assertions holds: 

1. The variational problems &25\) and V26}) are dual in the sense that = Doq. 

2. There exists a globally Lipschitz, anti- symmetric, convex-concave Hamiltonian H on R N x R N , such 
that the minimum in V25\) is attained at Lg, where H is the restriction of H to x O. 

3. There exists a measure preserving involution S such that the maximum in \26\) is attained. 

4- Moreover, for each H satisfying (2), there exists S satisfying (3) such that the following equation holds: 

u(x) =V 1 H(S(x),x) a.e. x E il. (27) 

2.1 Self-dual transformations 

We first introduce the following notion. 

Definition 2.2 Let fi be an open bounded subset of R w , and say that a measurable point transformation 
S : Cl — > O is self-dual if for every H £ L 1 (f2 x fi) such that H(x, y) = —H(y, x) for almost all (x, y) £ O x f2, 
we have 

/ H(x,S{x))dx = 0. (28) 
Jn 

In order to characterize these maps, we recall the following basic notion. 

Definition 2.3 A map S : Cl — > ft is said to be a measure preserving transformation if for every / £ L 1 (f2), 
foS£ L\ri) and 

/ foS(x)dx= / f{x)dx. (29) 
Jn Jn 
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By considering functions of the form H(x, y) — f(x) — f(y) where / £ i 1 (r2), one can easily see that self-dual 
transformations are necessarily measure preserving. 

The converse is however not true. Indeed, the map S(xi, X2) — (x2, —Xi) on R 2 is clearly measure preserv- 
ing. On the other hand, by considering the symplectic matrix J(x\,X2) = (—X2,xx) on R 2 , the function 
Hj(x,y) = (Jx,y) is clearly anti-symmetric. We therefore have that 

/ Hj(x,Sx)dx= / (Jx,Sx)dx=— / (Jx,Jx)dx = — / |a;|| 2 <ir. =/= 0. 
Jb Jb Jb Jb 

which means that S is not a self-dual transformation. 

More generally, recall that the linear and discrete analogue of measure preserving maps are the unitary 
matrices, i.e., those that satisfy UU* — U*U — I. If now U is a self-dual transformation, then one can easily 
see that U* = U and U is therefore a symmetric involution. This turned out to be true for more general 
transformations. 

Proposition 2.1 A measurable map S is a self-dual point transformation on fi if and only if it is both 
measure preserving and an involution, i.e., S 2 = I a.e., where I is the identity map on £1. 

Proof. Assuming that S is measure preserving such that S 2 — I a.e, then for every anti-symmetric H in 
L X {£1 x fi), we have 

H(x,S(x))dx = [ H(S{x),S 2 (x))dx = [ H(S(x), x) dx = - [ H(x,S{x))dx, 
; Jn Jn Jn 

hence J n H(x, S(x)) dx = 0. That is S is a self-dual transformation. 

Conversely, Let S be a self-dual transformation. We have seen above that it is necessarily measure preserving. 
Consider now the anti-symmetric functional H(x, y) — \S(x) — y \ — \S(y) — x\. We must have that 

= J H(x, S(x)) dx = J \S 2 {x) -x\dx, 

which clearly yields that S is an involution almost everywhere. 
An immediate application of Proposition 12.11 is that 

P > D 

Indeed, for any H 6 %, and any point transformation S on 17, we have 

L H (x, u(x)) + H(S(x), x) > (u(x),S(x)}. 
If now S is a measure preserving involution, we have that J„ H (x, S(x)) dx — 0, which means that 

Lh(x,u(x)) > / (u(x), S(x)) dx. 



2.2 Regularization of skew-symmetric functions 

Let f2 be a bounded domain contained in a ball Br centered at the origin with radius R > in K.^, and 
consider an arbitrary anti-symmetric function H 6 C(f2 x Cl). Our plan is to show that one can assign to H, 
a skew-symmetric convex- concave Lipschitz function H reg such that Lh < Lh on x Br. 
Note that for an elementary anti-symmetric function of the form H(x,y) = f{x) — f(y), where / is a 
continuous function on f2, one can easily show that H reg (x,y) = f**(x) — f**(y) where /** is a double 
conjugate of / defined as follows: 

f*{p) = sup{(x,p) - f(y)} and f**(x) = sup {(x,p) - f*(p)}. 
yen peB a 
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Note that these are not the usual Legendre transforms since the suprema are not taken over the whole space. 
The analogous regularization process for a general anti-symmetric function is more technical. 
We start by considering the class 

U- = {H G C(fi x n);H(x, y) < -H(y, x) for all x,y e tt}. (30) 

For each H G H- define Lh as in (|24|) and set 

= {L H ;Hen-}. 

For Lh G we define its (restricted) Fenchel dual L* H : R N x R N — > M by 

l h{i,v) = sup {(y,p) + - L H (x,p)}. 

x£Q,p£B R 

and similarly L*jf : R N x R N ->■ M by 

L H*(y,<i)= su p {(y,-p) + - l *h(p, x )}- 

For each convex function L : x R N ->• R, we shall define its B/f-Hamiltonian by 

H L (x,y) = sup {(x,p) - L(y,p)} 

p€B R 

Finally for each H G H, we define the convex-concave regularized -ff re9 of H by 

HL**(x,y) — Hl"(v,x) j. t 
H reg (x,y)= L » K l!J > 2 ^ , for all x, y G R N . 

We list some of the properties of H reg and Ln rscj ■ 
Proposition 2.2 Let H G T~L. The following statements hold: 

1. H reg is a skew-adjoint Hamiltonian on R N x R N whose restriction to Q x belong toU. 

2. Ln rag is convex and continuous in both variables and Ljj re < Lh on Q x Br. 

3. \L Hreg (x,p)\ < R\\x\\ +R\\p\\ + 5R 2 and \H reg {x,y)\ < R\\x\\ + R\\y\\ + 4R 2 for all x,y,p G R N . 
4- LH rrg and H reg are both Lipschitz continuous with Lipschitz constants less than ANR. 

To prove the above Proposition, we shall need the following two lemmas. 
Lemma 2.4 If H eH-, then hen iJ L « G H— 
Proof. For x, y G 1* we have 

H L' H *(x,y) = sup {(pi,x) - L*H{y,p l )}. 

Pl€ B a 

It follows from the definition of L*^ that 



l *h(v^Pi) = SU P \{zx,Pi) + {pi,y)-L H (jp2,Zi)\ 

\{zi,Pl) + (P2,y) - Sup {(Z 2 ,p 2 } + (P3,Zi) - L H (Z2,P3)\ 



sup 

Zi£Q,p2&B R k 22SO,p 3 gB H 



Sup inf I {zi,Pl) + (P2,V) - (Z2,P2) - (p3,zi) 

z 1 eQ,p 2 &B R z 2en,p 3 eB R >- 

+ sup {(z 3 ,p 3 ) - H(z 3 , z 2 )} \ 
z 3 en ' 

sup inf sup \ (zi,pi) + (p2,y) - (z 2 ,p 2 ) - (p3,zi) + (23, p 3 ) - H(z 3 ,z 2 ) ]■ 



(31) 
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Therefore, 

H L y{x,y) = sup _inf sup inf_ \ (pi, x) - (zi,pi) - (p 2 ,y) + (z2,p 2 ) + (p3,zi) - (z 3 ,p 3 ) 

pieB a z 1 en,p 2 eB R Z2 ^n, P3 GB R z 3£° 

+H{z 3 ,z 2 )). (32) 

Taking into account that H(z 3 , z 2 ) < —H(z 2 ,z 3 ) together with the fact that 

sup inf {...} < inf sup{...} 

yield 

H L* H »(x,y) < sup inf sup M_ sup I {pi,x) - (zi,pi) - (p 2 ,y) + (z 2 ,p 2 ) + (P3,zi) - (z 3 ,p 3 ) 
PieB a zien,p 2 £B Rp3 eB R z 3 eii Z2en ^ 

-H(z 2 ,z 3 )} 

= sup inf sup inf_ \{pi,x) - (zi,pi) - (p 2 ,y) + (p 3 ,zi) ~ (z 3 ,p 3 ) 

Pi£B a z 1 £(l,p 2 £B R P3 £B R z 3 eQ <- 

sup{(z 2 ,p 2 ) - H(z 2 ,z 3 )}\. 
z 2 en ' 

= sup inf sup inf _\{pi,x) - (z 1 ,p 1 ) - (p 2 ,y) + {p 3 ,z x ) - (z 3 ,p 3 ) + L„(z 3 ,p 2 )\ 

PiG-Br zi£(l,p 2 £B R P3 £B R z 3 eQ <- > 

Note also that inf Z3£ Q{...} < {...}| Z3=Zl , from which we obtain 

H L'*(x,y) < sup inf inf sup { (pi, x) - (zi,pi) - (p 2 , y) + (p 3 , zi) - (zi,p 3 ) + L H {zi,p 2 ) \ 
pieB R P2^B R z 1 en P3& B R L J 

= sup inf inf sup \ (pi, x) - (zi,pi) - (p 2 , y) + L H (z 1 ,p 2 ) \ . 
PieB R P2tB R z 1 en P3 eB R >- > 

We can now drop the term sup p3(EBj? since there is no p 3 in the expression in the bracket. Thus, 
H L **(x,y) < sup inf inf_ { (pi, x) - (21, pi) - (p 2 , y) + L H {z±,p 2 ) } 

pi<=B R P2^B R ZlG o I J 

= sup \(pi,x)+ inf mi{-(z 1 ,p 1 )-(p 2 ,y)+L H (z 1 ,p 2 )}\ 
Pl eB R >■ p 2 eB R zie o j 

= sup \(pi,x) - L* H (pi,y)\. 

pieB R > 

It follows from the definition of L*^ that L* H (p 1 ,y) + L* I ^(z,p) > (pi,z) + (y,p) for all z,p e CtxB R . Therefore 

H L * H *(x,y) < sup \(pi,x) - L* H ( Pl ,y)} 
PieB a *■ > 

< sup \{pi,x)- sup + (y,p) - L*g(z,p)}\ 



PiG-Bn *■ P eB R ,zeS2 

= sup inf _{(pi,x) - (pi,z) - (y,p) + L^*(z,p)} 
Pi&B R peB R ,zeii i. J 

< inf sup {(pi.z) - (pi,^) - (y,p) +L%(z,p)\ 
peB R ,zeQ Pl eB R ^ > 

< inf sup { (pi , x) - (pi ,x)-{y,p) + L*^ (x,p)\ 
peB RpieBa I J 

= mf R {-(y,p)+L* H *(x,p)}} 

= -H L , H ,(y,x). 

This proves that Hl** (x, y) < —Hl** (y, x) for all x, y G Cl. □ 



Lemma 2.5 If H £ H-. then L* H (p, x) < Lh{x,p) and L*^(x,p) < Ljj(x,p) for all x £ ft,p £ Br. 
Proof. For every (x,p) £ Cl x Br, we have 

L h(p, x ) = sup {(V,P) + (x,q) - L H (y,q)} 

= sup {(y,p) + (x,q) -sup{(q,z) - H(z,y)}} 

yeQ,q&B R 2GO 

= sup inf {(y,p) + (x,q) - (q,z) + H(z,y)}. (33) 

y en^eB R z ^ n 

It then follows that 

L* H (p,x) < sup inf {(y,p) + (x,q) - (q,z) - H(y,z)} 
y en, q eB R z ^ 

< inf sup {(y,p) + (x,q) - (q, z) - H(y,z)} 

z^ ye n,q£B R 

= inf sup{(y,p) + R\\x - z\\ - H(y, z)} 

z £ Q y£n 

= inf {R\\x — z\\ + Lh(z,p)} 

< L H (x,p). 

Thus, L* H (p, x) < Ljj(x,p). We now prove L*ff < Lh on Cl x Br. It follows from equality (|3"Tj) that for every 
{y,Pi) £ Cl x B R , 

L H{y,p-L)= sup inf sup I (zi,px) + {p 2 ,y} - {z 2 ,p2) ~ (P3,Zi) + (z 3 ,p3) ~ H(z 3 ,z 2 )\. 

Thus, 

L H(y>Pi) < -inf sup sup + (P2,y) - (^2,P2) - (P3,zi) + (^3,P3) -if(z3,Z2) I 

Therefore by considering Z2 = y and P3 = pi we get 

L *H*(y,Pi) < sup{(z 3 ,pi) - H(z 3 ,y)} = L H (y,Pi). 

This completes the proof. □ 
We now recall the following result from [5]. 

Proposition 2.3 let D be an open set in R d such that D C Br where Br is ball with radious R centered 
at the origin in R d . Let / : K d — > K and define f : R d -> R by 

f(y) = sup{(y,z) - f(z)}. 

zeD 

If f £ L°°(D), then f is a convex, Lip schitz function and 

\f(yi)-f(y2)\<dR\\ yi -y 2 \\, 

for all 2/1,2/2 £ R d - 

Lemma 2.6 If H £ then the following statements hold: 

1. \Ug(x,p)\ < R\\x\\ +R\\p\\ +3i? 2 and \H L %(x,y)\ < R\\x\\ + R\\y\\ + AR 2 for all x,y,p £ R N . 

2. L*jf and Hi,** are Lipschitz with Lip(H j,**) , Lip(L^) < 4NR. 
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Proof. It follows from the definition of Lh that Ln(y, q) > (y, q) on ft x Br. This together with 51 C Br 
imply that 

L *h(p, x ) = sup {{x,q) + (p,y) - L H (y,q)} 

< sup {{x, q) + (p, y) - (y, q)} 

< R\\x\\ + R\\p\\ + R 2 

for all x,p G l w . With a similar argument we obtain L*ff(x,p) < R\\x\\ + R\\p\\ + R 2 . We also have 
L h( x >p) = sup {(x,q) + (p,y) - L* H (y,q)} 

> (x, q) + (p, y) - L* H (y, q), ( for all [y, q) G Cl x B R ) 

> -R\\x\\-R\\p\\-R\\y\\-R\\q\\-R 2 

> -R\\x\\ -R\\p\\ -3i? 2 . 

Therefore < R\\x\\ + R\\p\\ + 3R 2 . The estimate for Hl*» can be easily deduced from its definition 

together with the estimate on L*^ . This completes the proof of part (1). 

For (2) set D = O x Br, then D C B 2 r where B 2 r is a ball with radius 2R in M 2Ar . Now assuming / = L* H 
in Proposition O we have that / = Ug. Therefore L*fi is Lipschitz in R 2N with Lip(L*^) < ANR. To 
prove that is Lipschitz we first fix y G and we define f y : M. N — >• R by / y (p) = L*^(y,p). Assuming 

D = Br G in Proposition 12.31 we obtain that the map x — > fy(%) — HL**{x,y) is Lipschitz and 

\H L »( Xl ,y) - H L »(x 2 ,y)\ <NR\\ Xl -x 2 \\ (34) 

for all x%,x 2 G M. N . Noticing that the Lipschitz constant NR is independent of y, the above inequality holds 
for all x\,X2,y G R . To prove HL^(x,y) is Lipschitz with respect to the second variable, let r > and 
2/1,2/2 G R w . Let pi and ^2 be such that 

(x,p 2 ) - L* H *(y 1 ,p 2 ) < H L **(x,yi) < (x,pt) +r, 

and 

{x,pt) - L*^(y 2 ,pi) < H L * H *(x,y 2 ) < (x,p 2 ) - L*^(y 2 ,p 2 ) + r, 

It follows that 

L*H(y 2 , P2 ) - L%( yi ,p 2 ) - r < H L ~(x,yi) - H L »(x t y 2 ) < L%(y 2 , Pl ) - L%( Vl , Pl ) + r 
L*^ is Lipschitz, thus, 

-4NR\\vx - y 2 \\ -r < H L ^{x,y x ) - H L ^{x,y 2 ) < WR\\ Vl -y 2 \\+ r 
Since r > is arbitrary we obtain 

-WR\\yi - 2/2 II < H L »(x,yx) - H L * H ,(x,y 2 ) < WR\\ yi -y 2 \\. 
This together with ( [34]) prove that H^*^ is Lipschitz and Lip(Hl**) < ANR. □ 

Proof of Proposition I2T21 It is easily seen that H reg (x,y) — —H reg (y,x) for all x, y G R . 

Now note that by definition 

HL**{x,y) = sup - L*x(p,y)}, 

pes R 

and therefore for all ?/ G R w , the function x —> H^** {x, y) is convex. We shall show that for all x G M. N the 
function y — > (x, y) is concave. In fact we need to show that 

y -> -H L **(x,y) = inf {L*^(p,y) - (x,p)} 
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is convex. For that consider A £ (0,1) and elements 2/1,2/2 G 1*. For every a > —Hi J ^(x,yi) and 6 > 
^H^^(x,y 2 ), find pi,P2 G such that 

-H L , i *{x,y 1 ) < L*H(p 1 ,y 1 ) - (x,px) < a and - H L **(x, y 2 ) < L*£ (p 2 ,y 2 ) - (a:,p2> < b. 

Now use the convexity of the ball Br and the convexity of the function L*^ in both variables to write 

-H L » H *(x,\ yi + (1 - X)y 2 ) = inf {L* H *(p, \ Vl + (1 - A)y 2 ) - (x,p)} 

< L* H *(\ Pl + (1 - A)p 2 , At/i + (1 - X)y 2 ) - {x,\ Pl + (1 - X) P2 )} 

< X(L*H(pi,yi) ~ (x,pi)) + (1 - A)(i^*(p 2 ,2/2) - (£,£2))} 

< Aa + (1-A)6. 

This establishes the of concavity of y — > H^{x 1 y). It then follows that H reg (.,y) is convex and H reg (x, .) 
is concave on l w , which completes the proof of part (1). 
We now prove part (2). Let (x,p) G x Br. We have 

L Hrcg (x,p) = snp{(y,p) - H reg (y,x)} 
yen 

H L **(y,x)-H L *.(x,y) 

= sup{(2/,p) s ^ 5 } 

yen z 

< sup{(y,p) -H L * H *(y,x)} (by LemmaED]^" (2/, x) < -H L % [x, y)) 
yen 

= sup {(y,p) - sup {(?/,g} - L*h (x,q)}} 
yen qeB R 

= sup inf {{y,p)~ (y,q) + L%{x,q)} 
yen q£B R 

< inf sup { {y, p) - (y, q) + (1, g)} 

< sup {(y,p) - (y,p) + L*h (x,p)} 
yen 

= L* H *(x,p). 

Thus, Ln rc (x,p) < L*^{x,p). It also follows from Lemma 12.51 that L*^(x,p) < Lh(x,p) from which we 
obtain Ln rcg < Lh on i7 x Br. 

Hl** { x :V)~ Hl** {y,x) 

By the dchnition H reg (x, y) — — ^ . Therefore, by Lemma lzToI H reg is Lipschitz and Lip(H reg )) = 

Mp(H L ~) < ANR and also 

\H reg (x, y)\ < 2 S < R\\x\\ + R\\y\\ + 4R 2 . 



The corresponding results for Lh TS9 follow by the same arguments as in parts (2) of Lemma 12.61 and the 
bound for H reg . □ 



2.3 Proof of Theorem 1231 

We first show that the minimization problem (f2"5t has a solution. Let Br be a ball such that f2, C Br. 
Let {#"} be a sequence in H such that Ln n is a minimizing sequence for P M . It follows from Proposition 
12.21 that L/fn < L# n on x i?R and therefore Ljjn is still a minimizing for P^. It also follows from 
Proposition I2.°2I that £ff» e and -ff™ eff are uniformly Lipschitz with Lip(H^ eg ) , Lip(L ji» ) < 4iVP and also 

I H?e g (x,y) I < P|M| +P||2/|| +4P 2 and |L ff?cg (z,p)| < P||.x|j + P||p|| + 5i? 2 
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for all x,y,p G Mr. By Arzela-Ascoli's theorem, there exists two Lipschitz functions H, L : R N x R N — > K 
such that H" eg converges to H and Ln^ Bg converges to L uniformly on every compact set of M. N x M. N . This 
implies that H G %. Note that 

LH ?eg (x,p)+H? eg (y,x) > (y,p), 
for all x,p 6 R N and y G from which we have 

£(z,2>) > (y,p) - H(y,x), 

for all x,p £ M. N and y G f2. It implies that L# < L. Let -ff res be the regularization of H defined in the 
previous section. Set — H reg and = Lh^ . It follows from Proposition 12.21 that Ljj^ < Ljj on 
Cl x Br, from which we have 



■P-y 



/ Lh{x, u(x)) dx — I L oc (x, u(x)) dx. 
Jn Jn 



□ 

For the rest of the proof, we shall need the following two technical lemmas. 

Lemma 2.7 Let Hoo be the Hamiltonian obtained above. For each x G Cl, define f x : WL N — > K by f x {y) = 
Hoo(y,x). We also define f x : l w — > K U {+00} by f x (y) — fx{y) whenever y G and +00 otherwise. Let 
{fx)* be the standard Fenchel dual of f x on M. N , in such a way that (f x )*** = (f x )* on Wi N . Then we have, 

1- fx = {fx)** = fx on tt, and 

2. Loo(x,p) = sup{(z,p) - {fx)**{z)} = sup {(z,p) - {f x )**{z)}. 

2ES1 zGR" 

Proof. (1) Since {f x )** is the largest convex function below f x we have and f x < {f x )** < f x , from which 
we obtain f x — {f x )** — f x on fl. 
For (2), we first deduce from (1) that 

Cfx)*{y) = Cfx)***{y) = sup{( 2 ,y)-(/ x )**(z)} 



> sup {(«,»> -(/»)**(«)} 

> sup{(z,y> -(/*)**(*)} 

= su.p{(z,y) - f x (z)} 
zen 

= sup{(z,y) - f x {z)} 
zen 

= {fx)*{y), 

from which we have the desired result. □ 
Lemma 2.8 Let H G H. For each r£l and A > define 

L r>x (x,p) := S up{{z,p)-{f x r{z)-X^f + X^f-rH{z,x)}, 

L x {x,p) := S np{(z,p)-{f x r{z)-X^f- + X^f}, 
L r {x,p) := sup{(z,p) - H x (z,x) - rH(z,x)}. 

Then the following assertions hold: 

1. For every {x,p) G M. N x Mr , we have lim L\{x,p) — Loo{x,p) and lim L r \{x,p) — L r {x,p). 
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2. For all x E R N , the function p — > L x (x,p) is differentiable. 

3. For every (x,p) ER N xR N , we also have lim L r.^,p)-Lx(^,p) = H(\7 2 L x {x,p),x). 

'I 7-0 r 

Proof. Yosicla's regularization of convex functions and (1) of Lemma l2~7l yield that 

lim L rt \(x,p) = sup{(z,p) - (f x )**(z) - rH(z,x)} = sup{{z,p) - H^z^x) - rH(z,x)} = L r (x,p). 



We also have 

lim L x (x,p) = sup {(z,p) - (f x )**(z)}, 

which, together with (2) of Lemma [2~71 yield limA-j-o L x {x,p) — L^x^p). 

(2) follows from the fact that the Yosida regularization of convex functions are differentiable. 

(3) We let z rA E fi and z' r X € R N be such that 

L r ,x(x,p) < (z r ^p)~(f x )**(z rA )-X^^ + X^-rH(z rA ,x)+r 2 , 

Therefore, 

- H(z> x ,x)-r< L rA^P)~L x (x : p) ^ _ + r 

r 

Note that by the definition of L x , we have sup r£ [_ 1 ji \\z' r x \\ < oo. Suppose now that, up to a subsequence, 
z r,\ z x E Cl and z' r A — > z' x as r — > 0. This together with the definition of L r x and L x imply that 

L x (x,p) = (z x , P ) (f x r(z x ) - X l ^t + X M1 = {z > X}P) - A^ + A^, (36) 

from which we obtain 

z\ = z' x = V 2 L x (x,p) € fi. (37) 

Therefore, it follows from (|35j) that 

L rX (x,p) - L x (x,p) 
km = H(V 2 L x (x,p),x). 

i — >o r 

□ 

End of the proof of Theorem l2.lt For each A > 0, x E Cl and p E M. N , we define S x (x,p) = \72L x (x,p). It 
is easy to see that S x (x,p) — > So(x,p) where So(x,p) is the unique element with minimal norm in d2L oa (x,p) 
(see Proposition 1.3 in pQ). Set S(x) — Sq(x, u(x)). For each r > 0, A E [0, 1] and x E fi, define 

L r x (x,u(x)) - L x {x,u{x)) 

rj r {\,x) = . 

r 

Note that the function r — >• L r X (x, u{x)) is a convex function because it is supremum of a family of linear 
functions. Thus, for fixed (a;, A) E £1 x [0,1], the function r — > rj r {\,x) is non-decreasing. Setting r]o(X,x) 
to be H(S x (x),x) for A > and 770(0,2;) = H(S(x),x), we have that both functions A — > r/ r (A,x) and 
A — ^ 770 (A, x) are continuous. It follows from Dini's Theorem, that for a fixed x, ?7 r (A, x) converges uniformly 
to r]o(X, x) as r — > with respect to A E [0, 1]. 

Note also that thanks to (|37p we have that S x , S : Ct — > fi. We shall now show that J n H(S(x),x) dx = for 
every H E %, meaning that S is indeed a self-dual point transformation. Indeed, by Fatou's lemma we have 

lim / H(S x (x), x) dx = / H(S(x),x)dx. 
A->0./o In 
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It follows from (|55"|) that 

It follows that 

H(S(x),x) dx 



L r< x(x,p) - L x (x,p) 



< \\ H \\l°-(B r xB r ) + 



lim lim 

q A^O r^0+ 



L r< x(x,u(x)) - Lx{x,u{x)) 



dx 



I lim lim ri r (X,x) dx 



lim lim r? r (A, x) dx 

Q r^0+ A^O 

lim r/ r (0, x) dx 
n r->0+ 



lim 

r->0+ 



77 r (0, a;) dec 



(due to the uniform convergence) 



(due to the monotone convergence theorem) 



lim 

r-s-0+ 



L r (x,u(x)) - £oo(a:, u(x)) 



dx 



> 0, 



from which we have H(S(x), x) dx > 0. 

By the same argument considering r —> 0~, one has J„ H(S(x),x) dx < and therefore the latter is indeed 
zero as desired. 

It follows from the fact that S(x) £ d2L 00 (x,u(x)) together with (f x )** being the Fenchel dual of L with 
respect to the second variable (in view of Lemma |2~T)) that u(x) £ d(f x )** (S(x)). By considering Theorem 
3.11 in the Appendix, assume that ft' is a dense subset of Br such that for each z e Of, \7iH 00 (z,x) exists 
for all x € Vt. Define Slo = 5 _1 (r2' \ dil). Since S is measure preserving we have that fio is dense in Cl. We 
also have for each x G Clo, that ViHoq(S(x),x) exists. Since (/»)**(•) = Hoc{-,x) on we obtain 



w(x) = ViH oc (S(x),x), 



for all a; £ SIq- 



To complete the proof of Theorem l2.ll it remains to show that Poo = -Doo- We already know that Poo > D c 
To prove the equality it suffices to notice the following: 



p — 



L 00 (x, u(x)) dx 



L rx (x, u{x)) dx + / H 00 (S(x), x) dx 
Jn 

Loa(x,u(x)) dx + / (f x )**(S(x)) dx 



(u(x), S(x)) dx < D c 



2.4 Remarks on the uniqueness of the decomposition 

We have seen in example (14) that one cannot expect uniqueness of the involution S in the above decom- 
position of a given vector field u. We now complete the part of the proof of Theorem 11.11 which gives the 
uniqueness of the involution. (|SJ) on u. 

Proof of Theorem 11.11 Assume first that u(x) = ViHi(SiX,x), for some Hamiltonian Hi and some 
selfdual transformation Si. We shall show that (Hi, Si) is an "extremal pair" (i.e., where Poo and 6*00 
are attained), and that u(x) = ViPoo(<Si£, %), where Poo is the optimal Hamiltonian constructed above. 
Indeed, let L be the Fenchel-Legendre dual of Pi with respect to the first variable. We have that Lh 1 < L 
on R N x fi. It follows that 

(u(x), Si(x)) < Lh ± (Si(x),u(x)) + Hi(Six, x) < L(x, u(x)) + Hi(Six, x) = (u(x), S(x)) 
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from which we have 

(u(x),Si(x)) = L Hl (x,u(x)) + Hi(Six,x), 

and 

/ (u(x), Si(x)) dx = I Lh 1 (x, u(x)) dx. 
Jn Jn 

On the other hand we have 

/ (u(x), Si(x)) dx < Doc — Poo < / L Hl (x,u(x))dx, 
Jn Jn 

which yields 

/ (u(x), Si(x)) dx = Doc — Poo = / L Hl (x,u(x))dx. 
Jn Jn 

Now we can show that u{x) — Vi-ff 00 (5 , ix, x). In fact, 



(u(x), Si(x)) dx = / Lh 1 (x,u(x)) dx — Poo — I L oa (x, u(x)) dx + / Hoo(Si(x), x) dx, 
n Jn Jn Jn 

which implies that 

(u(x), Si(x)) = Loo(x,u(x)) + Hoo(Si(x),x) 
a.e. on f2, and hence the desired result. 

Assume now that the function x — > (u(x),yi — j/2) + i u (yi) — u (y2),x} has no critical point unless when 
Vi =2/2- Suppose 51,5*2 are two transformations such that for i = 1,2, we have 

(u(x), -u(Six)) — ViHi(SiX,x). (38) 

We shall show that S\ = S2 a.e. on fi. Note first that the previous argument gives that 

(u(x), -u(Six)) — ViHoo{S t x,x). (39) 

Note also that the function x — > L x (x,u{x)) is locally Lipchitz and therefore is differentiable on a subset 
VLq of full measure. We now show that 5*1 = 5*2 on fio- 
Indeed, for any x 6 fio, h = is a minimum for the function 

h — > Loo(x + h, u(x + h)) + Hoa(Si(x), x + h) — (u(x + h), Si(x)). 

This implies that 

V 2 i?oo(S , i(a;), x) - Du{x)Si(x) = --^rL^x + h, u(x + h)) h=0 = V 2 H o('5'2(a;), x) - Du(x)S 2 (x), 

dh 

from which it follows that 

Du(x)(S 2 (x) - Sx(x)) = u(5i(a;)) - u(S 2 (x)). 
The hypothesis then implies that S\(x) = S2 (x), and S is therefore unique. 

Remark 2.9 Note that the function 

h — > Loo(x + h, u(x)) + H tx (S(x),x + h) — (u(x), S(x)). 
has a minimum at h = 0, from which we have 

u(S(x)) = -V 2 i?oo(5(x),x) G diLoo(x,u(x)). 
Since, on the other hand, we have S(x) S d2L 00 {x,u{x)), one obtains 

(u(S{x)),S(x)) edLoo(x,u(x)). 
Now under the hypothesis (0) that ensures uniqueness, the above inclusion becomes 

(u(S (x)) , S (x)) = VLoo(x,m(x)) a.e. x G f2. 
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Suppose now that u is a monotone map. We shall show that it satisfies condition ©, This will then yield 
(3) of Theorem II .11 since u is then a.e., differentiable and the decomposition holds with S being the identity, 
according to the theorem of Krause. 

Indeed, any critical point x of the function x — > (u(x),yi — y 2 ) + (u(yi) — u(y 2 ),x) satisfies 

Du(x)(yi - y 2 ) + u(y\) - u(y 2 ) = 0, 

hence 

(Du(x)(yi - y 2 ), yi - y 2 ) + (u(y%) - u(y 2 ), y\ - y 2 ) = 0. 

Since both terms are non-negative, they are equal to zero. If u is strictly monotone, this cannot happen 
unless yi = y 2 . 



3 Appendix 

Theorem 3.1 Let H be a skew- symmetric finite convex-concave function on R N x M. N such that for some 
A > 0, it satisfies 

\H(x u yx) - H(x 2 ,y 2 )\ < A[[i x -x 2 \\ + A\\yi - y 2 \\ for all (x t , yi ), (x 2 ,y 2 ) eR w x R N . (40) 

Let A C M. N be a closed ball and let B C be a compact subset with non-empty interior. Then, there exists 
a dense subset A' of A such that for each x 6 A', ViH(x,y) exists for all y £ B. 

This result is actually a particular case of a more general result established in [5] , where the same conclusion 
is established for finite convex-concave functions on M" x W n with n ^ m and without condition (|4T)]) . For 
the special case n = m = N considered in Theorem 13.11 the proof can be shortened and we shall provide 
here a sketch for the reader's convenience. We shall need a few preliminary results. 
The following definition and theorem can be found in [2]. 

Definition 3.2 A sequence {f n } of (scalar-valued) functions on an arbitrary set X is said to converge to 
f quasi-uniformly on X, if {/„} converges pointwise to f and if for every e > and L € N, there exists a 
finite number of indices ni,n 2 , > L, such that for each x G X , at least one of the following inequalities 
holds: 

\f ni (x) ~f(x)\ < e, i = l,2,...,k. 

Theorem 3.3 If a sequence of functions on a topological space X converges to a continuous limit, then the 
convergence is quasi-uniform on every compact subset of X. Conversely, if the sequence converges quasi- 
uniformly on a subset of X, then the limit is continuous on this subset. 

For (x, y) £ ffi.^ x R w , the one sided directional derivative of H at (x, y) with respect to (u, v) is defined as 
the limit 

H(x + Xu,y + Xv) - H(x,y) 
VH (x, y)(u, v) = hm 

A^0+ A 

provided such a limit exists. It is standard that the directional derivatives 

rji u „ x ,• H(x + Xu,y)-H(x,y) 

ViH(x,y)(u) = hm 

and 

„, w x r H(x,y + Xv) - H(x,y) 

V 2 H(x,y)(v) = hm 

exist. The following result is due to T. Rockafellar |10) . 

Theorem 3.4 Let H be a convex- concave function on R N x WL N . Let C x D be an open convex set on which 
H is finite. Then for each (x,y) G C x D, V ' H(x,y)(u,v) exists and is a finite positively homogeneous 
convex-concave function of (u, v) on R N x R N . In fact, 

VH(x, y)(u, v) = V x H(x, y)(u) + V 2 H(x, y)(v). 
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For each A > 0, define the following functions on R N x R N . 

H(x + Xu, y + Xv) — H(x, y) 



H\(u,v) 
H\(u,v) = 



X 

H(x + Am, y + Xv) — H(x, y + Xv) 



X 

Hl(u) = H\(u,0) and H^(v) = H\(Q,v). Note that H\ and H\ are monotone, so that by Dini's Theorem, 
both H\{u) and H\(y) converge uniformly on compact subsets of Ht N to V\H(x,y)(u) and V%H(x,y){v) 
respectively. We have the following properties for H\. 

Proposition 3.1 The following statements hold: 

1. H\(u, v) converges uniformly to ViH(x, y)(u) + \?2H(x, y)(v) on compact subsets Ax B ofM. N x Mr. 

2. If for some u £ R N we have X/\H(x, y)u = —ViH(x, y)(—u) then for each v E B. 

lim A ^ + Vi-ff (x, y + Xv)u = V\H{x 1 y)u uniformly on B. 

3. H\(u,v) converges uniformly to ViH(x,y)(u) on compact subsets Ax B ofM. N x R*. 
Proof. We first show that for each e > 0, there exists Ao > 0, such that for all < A < Ao, we have 

H x {u,v) < \7 1 H(x,y)(u) + \7 2 H{x 1 y){v) + e, for all (u,v) G Ax B. 

Then by a dual argument we have 

H x {u,v) > ViH(x,y)(u) + V 2 H(x,y)(v) - e, for all (u,v) e Ax B, 

from which we obtain the desired result in part (1). The difference quotient in the function H\ can be 
expressed as 

H (x, y + Xv) — H(x, y) H(x + Xu, y + Xv) — H(x, y + Xv) 
A + A ' 

where the first quotient converges uniformly to V2-ff(a;, y)(v) on B. Since H^(u) converges uniformly to 
V ' \H{x,y){u) on A, there exists a > such that 

H(x + au,y) - H(x,y) 

< ViH(x, y){u) + e. 

a 

Since H is Lipschitz with Lipschitz constant A > 0, for every uefi, we have 

H(x + au,y + \v)-H(x,y + Xv) < y ^ , + £ + 2XA\\v\\ ^ 
a ''2a 

Let Ao be small enough such that 2A °^^ < e/2 for all v 6 B. For each < A < min{Ao, a} we have 

H(x + au, y + Xv) — H(x, y + Xv) 



ViH( Xl y)(u) + e > 



> H(x + Xu,y + Xv) - H(x, y + Xv) 



from which part (1) follows. 

We know prove part (2). Note first that 



H(x + Xu,y + Xv) - H(x, y + Xv) 

^ > Vi-H {x, y + Xv)u, 

from which together with (|4ip . we have 

V 1 H(x,y)(u) + e>W 1 H(x,y + Xv)u, (42) 
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for every < A < min{Ao, a} and v G B. By a similar argument we have 

- Vi#(x, y)(-u) - e < -ViH{x, y + Xv)(-u), (43) 
It follows from (gSJ) and (@3J that 

-Vii/(x,y)(-u) - e < -V t H(x,y + \v)(-u) < ViH(x,y + Xv)u < V 1 H(x,y)(u) + e, 
and the result follows due to assumption V\H{x,y)u = —ViH(x,y)(—u), 

Part (3) follows from the fact that H\(u,v) — H\[u,v) — H^(v), □ 

Proposition 3.2 Fix u G R. . There exists A' C A with A' dense in A, such that for each x G A' , there 
exists a dense subset of B, say B x ^ u , such that for all y G B xu we have X? \H(x,y)u = —\/iH(x,y)(—u). 

Proof. Define F : A — > K by F(x) = J B H(x, z) dz. Note that F is convex and therefore there exists a dense 
subset A' C A on which F is differentiable. For every x G A' we have 

f „ [If (a: + Au, z) - H(x, z)] dz j [H(x + Xu,z) - H(x,z)] dz 
b (x)u — hm — £ - J = i = hm — £ - J = i . 

A-V0+ A a^o- A 

due to Lebesgue monotone convergence theorem we have 

. f H(x + Xu, z) — H(x, z) f H(x + Xu,z) — H(x,z) 

b [x)u = / hm dz = hm dz, 

J B \^o+ A Jb^o- X 

from which we have 

/ [V 1 H(x,z)u + V 1 H(x,z)(-u)]dz = 0. 

JB 

Since the integrand is nonnegative there exists a dense subset B XtU of B such that 

ViH(x, z)u + V\H(x, z)(— u) = 0, for all z G B XtU . 

□ 

Proof of Theorem 13.11 Let A' be as in the above Proposition. Fix x G A'. We shall show that for all 
u G R N and y G B, we have ViH(x,y)u + ViH(x,y)(—u) = from which we obtain ViH(x,y) exists for 
all y G B. Fix u G M w and define f(y) = V \H(x,y)u. We first show that / is continuous on B. Note first 
that f(y) = lim„^ 00 f n (y) where 

H(x + X n u,y)-H(x,y) 
fn(y) = : , 

and A„ = 1/n. We shall show that /„ converges quasi-uniformly to / on B. Fix e > and L G N. It follows 
from Proposition 13.21 that there exists a dense subset B x . u of B such that 



Vifl"(a:, y)u + VxHix, y){-u) = 0, for all y G B x>u . 
For each y G B X}U , it follows from Proposition 13. II that there exists n y > L such that 



H(x + X ny u,y + X„ y v) - H(x,y + X Uy v) 



X„ 



-M 



e 

<2' 



and 

|Viff(x,j/ + A nt/ u)u- /(y)| < -, 

for every v G -B. This implies that 

|/n v (y + A nyU )-/(T/ + A nv t;)| <e, (44) 
for all v G -B. Define ?7 y = {y + A ra «; u G B}. Since B KjU is dense in B we have 

B c Uy£B x , u int(U y ). 
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B is compact and therefore there exist y\,yi, ■■■,yk & B XfU such that B C D^ =1 int(U yi ). This together with 
(1441) implies that /„ converges to / quasi-uniformly on B and therefore / is continuous. 

Since / ic continuous and \J\H(x, y)u + VxH(x, y)(—u) = for almost all y £ B, we indeed have 

ViH(x, y)u + ViH(x, y){-u) = 0, 
for all y £ B. This completes the proof. □ 
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